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THE PROPER WEIGHTING FUNCTION 
FOR RETRIEVING TEMPERATURES 
FROM SATELLITE MEASURED RADIANCES 


One class of methods for converting satellite measured radiances into atmos- 
pheric temperature profiles, involves a linearization of the radiative transfer 
equation - e.g.. 


AR = 



( 1 ) 


where AT. is the deviation of the temperature in layer i from that of a refer- 
ence atmosphere, AR is the difference in the radiance at satellite altitude from 
the corresponding radiance for the reference atmosphere, and W, is the dis- 
crete (or vector) form of the T-weighting (i.e., temperature weighting) function 
W(P), where P is pressure. The top layer of the atmosphere corresponds to 
i = 1, the bottom layer to i = s - 1, and i = s refers to the surface. Lineariza- 
tion in temperature (or some function of temperature) is at the heart of all linear 
or matrix methods. The question we raise here is: What is the weighting func- 
tion that should be used in Eq. (1) ? 

Methods based upon statistical regression determine W empirically, but those 
methods that involve direct inversion (see Fleming and Smith, 1972, for an excel- 
lent review) use the expression 


W<*> 


dB(T) 

dT 


At. 

1 


where B(T) is the blackbody radiance at temperature T, A/. t._j - r., 

is the transmission from the top of the atmosphere to the bottom boundary of 
the ith layer for i < s, r 0 = 1, and r s = 0. All quantities are calculated for 
the reference atmosphere. The continuous form of the atmospheric part of 
W/ l) is 


dT 


T(P) 


dr(P) 
d 1 n P ’ 


where InP is the height parameter. This expression for the weighting function 
is correct if t(P) is independent of the temperature profile T(P). Otherwise, 
there is an additional term that comes from differentiating r with respect to 
temperature. (In performing the integration, one must note that r(P) is a func- 
tional of T(P') for 0 ^ P' ^ P, which, where necessary to be explicit, will be 
denoted by r [T(P')|g ; Pi). The proper linearization of the radiative transfer 
equation, when r depends upon the atmospheric temperature profile, leads to the 
proper T- weighting function 


w = w< n + w (2) 

1 11 

where 


x — i dAx. 

W< 2 > = ) B(T ) . 

1 Z_j J dT. 

_ . 1 
) * 1 

In continuous form, W(P) = W (1) (P) + W (2) (P), where 


( 2 ) 


W< 2 >(P) = + B(T S ) 


dpT [T(P;| 0 , P S 1 
[dT (P) d In P] F 



B(T (P* )) 


dpdr [T (P ) I Q * ; P‘1 

[dT (P) d In P] F In P' 


d In P', 


(3) 


where d F denotes the functional derivative. The first term contains the func- 
tional derivative of the transmission to the surface with respect to T(P); in the 
integral the transmission has been differentiated twice: first a (logarithmic) 
derivative with respect to P’ and then a functional derivative with respect to 
T(P). These expressions are derived in the Appendix. 

With the T-weighting function properly defined, its physical significance is 
readily apparent. W(P) is the change in the radiance measured at the top of the 
atmosphere per unit change of temperature over a unit extent of the height 
parameter In P. Thus, 


f p s 

AR = W S AT S + W(P) AT(P) d In P. 

Jo 

which is the continuous form of Eq. (1). Of course, this equation is strictly cor- 
rect only in the limit A T - 0 and is practically correct whenever the non-linear 
terms in the relation between temperature and radiance can be neglected. The 
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corresponding equation using W*(P) can be considered practically correct only 
with the additional limitation that the dependence of transmission on temperature 
can be neglected. 

A comparison of W (l) (P) with W(P) for typical sounding channels (NIMBUS 6 
HIRS channel 5, centered at 717 cm -1 , and channel 13, centered at 2244 cm 1 ) 
is shown in Fig. 1, calculated for the U.S. Standard Atmosphere (mid-latitude, 
spring/fall). The differences are significant. For the 717 cm -1 channel, the 
integral of W (n (P) from 0 to P s is 1.09 ergs/(cm • sec • str • °K), while the 
corresponding integral of W(P) is 0.79. Thus, the neglect of W* 2) (P) would 
overestimate by 38% the effect on the radiance of a uniform change in atmos- 
pheric temperature. For the 2244 cm -1 channel, the overestimate of the radi- 
ance due to neglect of W (2) is about half as much, 18%. For non-uniform 
temperature changes, the error can be larger or smaller. 

In addition to reducing the sensitivity of the channel to temperature changes, the 
additional term a^o shifts the position of the weighting function. This is illus- 
trated by W REN (dashed curve on left hand side of Fig. 1) which is W renormalized 
so that its peak value is the same as W M ) . The peak of W is shifted downward 
from W (1) by about 60 mb, it is narrower above the peak, and has a larger tail 
in the stratosphere. 

The effect of the W (2) term on temperature retrievals can best be assessed by 
considering W ( 1 ) to be the T-weighting function incorporated into the retrieval 
process and corresponding to a channel for which measurements are desired. 

The actual measurements, however, are different, corresponding to the T- 
weighting function W. That difference can be treated as a measurement error: 

s 

\R_ = r W (2) AT. . 

E L ■ 1 

i * 1 

The error thus depends upon how far the solution is from the initial guess. For 
a good initial guess, such as that obtainable from forecasts, a typical value of 
AT in the troposphere (below 200 mb) is ~ 2 to 3°K. Under those circumstances, 

L'.Ri = 2.5 

E 

which has the value .50 and .0029 ergs/(cm • sec • ster) for the 717 cm" 1 and 
2244 cm" 1 channels, respectively. At 243°K, the brightness temperature asso- 
ciated with these channels for the standard atmosphere, these radiance errors 
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correspond to brightness temperature errors of .43 C K for the 717 cm” 1 channel 
and .21®K for the 2244 cm" 1 channel. 

The above errors, treated as radiance errors, are larger than the instrument 
noise usually specified for these channels. Furthermore, these errors tend to 
be consistent over a synoptic scale grid (~500 km), so that averaging the meas- 
ured radiances within a grid will not reduce this error. Finally, the statistical, 
or minimum-RMS, method and the closely related minimum information method 
(reviewed by Fleming and Smith, 1972), which minimize the effects of radiance 
errors on the retrieved profile, assume that radiance errors are independent of 
atmospheric temperature, but that is certainly not the case with the error asso- 
ciated with the weighting function definition. Thus, with any of the linear methods - 
such as the statistical method, the minimum information method, as well as the 
Backus -Gilbert method (Conrath, 1972), which rely on calculated (as opposed to 
empirically determined) weighting functions — the neglect of the W ( 2 5 term in 
the weighting function may seriously affect the accuracy of the results. 

For weighting functions associated with very narrow spectral intervals, essen- 
tially monochromatic, the relative importance of W ( 2) can be even larger. An 
example of such a channel, centered at 2386.88 cm" 1 and 0.3 cm -1 wide, is shown 
in Fig. 2. The effect of W (2 \ which is usually negative for reasons explained 
below, is so strong that the proper T-weighting function is actually negative above 
400 mb. It simply means that increasing the temperature in the layers above 
400 mb will actually reduce the radiance in this channel. The "radiance error" 
due to the neglect of W (2) would be .014 ergs/(cm . sec • ster) for At <= 2.5°K. 

This corresponds to an error of approximately 1°K in brightness temperature at 
257 °K, the brightness temperature for this channel with a standard atmosphere. 

The factors which determine the sign and magnitude of W ( 2) (P) can be brought 
out by inserting into Eq. (3) the general form for transmission 


r(P) = r [T (P' ) | p ; P] = <t v [T(P')|P; P]> 


r„(P) 5 T v [T (P* ) I o ’ A = exp 


-f 

•>0 


K V (T(P'), P') d In P' 
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where K v — the volume absorption coefficient at wave number v, times the local 
scale height — is a point function of pressure P’ and the local temperature T(P’). 
(The angular brackets denote averaging over the spectral response of the channel.) 
Carrying out the differentiation, with care in taking functional derivatives, leads 
to the following expression 


W(2)( P ) 


“K v (T(P), P) 
dT(P) 


B(T(P))t„(P)-B(T s ) t v ( P s ) 


+ 



B(T(P' )) 


dr v (F) 

d In P' 



( 4 ) 


The first term in the square brackets corresponds to the diagonal of the matrix 
d /Sr./dT i appearing in Eq. (2), which expresses W (2) in vector form. Since 
d r^/dlnP is negative, it follows that the diagonal term is always of opposite 
sign to the off-diagonal elements P ? < P. (The off-diagonal elements on the 
other side, P’ < P, are of course zero.) Assuming no temperature discontinuity 
at the surface, Eq. (4) can be integrated by parts to yield 


W< 2 >(P) = - 


'dMT(P), P) 
dT(P) 


f P % ( p.) damped in F 
1 dlnP' / 


( 5 ) 


We are now in a position to explain the sign of W ( 2) (P). For the numerical 
examples considered here (Figs. 1 and 2), dK v /dT > 0 because of the location 
of these channels on the wings of the C0 2 absorption band; for the dominant 
absorption lines, the rotational quantum numbers of the lower states involved 
in the transitions are high and the rotational energies are > kT. It then follows 
that the occupation probabilities for those states increase with increasing 
temperature, resulting in an increase of absorption with increasing temperature. 
With dK„/dT > 0, the sign of W (2) (P) depends upon the temperature gradient. 
When dT(P)/dlnP > 0, as is generally the case in the troposphere, W (2) (P) 
becomes negative. 

The magnitude of W (2) (P) is very sensitive to the ratio of the rotational energy 
E to kT. For E > > kT (e.g., on the wings of the vibrational-rotational 
absorption bands) dK^/dT is large; but for E rot < kT (e.g., in the Q-branch) 
dK„ /dT is small in magnitude and could be positive or negative. Eq. (5) also 
reveals that W (2) (P) is approximately proportional to the temperature gradient, 
and for an isothermal atmosphere, it is zero. 
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We conclude, therefore, that the contribution of W (2) (P) to the proper weight- 
ing function depends upon the spectral position of the channel within the absorp- 
tion band and upon the temperature gradient. The largest effect is in the tropo- 
sphere, as in the examples shown here, and one would expect the effect to be 
smaller in stratospheric channels (in fact, it turns out to be negligible in the 
Q-branch channels). 

It is to be noted that the considerations in this note apply as well to methods 
which linearize the radiative transfer equation with respect to the Planck func- 
tion instead of temperature. In that case, the weighting function usually employed 
is dr(P)/d In P, but the proper weighting function should be 


v p ) = 


dT(P) 
d In P 


+ 


1 

dB dT 


W< 2 >(P); 


r 


L'l 


and the change in radiance is related to the change in Planck radiance ^B(P) at 
each level by 


AR = 



W H (P) AB(P) d In P. 


In more general applications — e.g., nonlinear retrieval methods in which the 
radiative transfer equation is iteratively applied to a temperature profile which 
undergoes correction at each iteration step — functional deri /atives can be used 
in a Taylor series expansion to relate the transmission function for one temper- 
ature profile T(P) to the transmission function for a corrected temperature 
profile T(P) +AT(P). For example, if the original transmission function is 
r (P) then the transmission for the corrected temperature profile is given by 


T(P) + 


f 

Jo 


1 r r _ 

2 Jo Jo l- 


d F T[T(P')|; ; P] 

— - - AT(P' ) d In P' 

[dT (P' ) d In P'] F 

( ip 7 1 T (P * ) I © ; P] 


1T(P" ) d In P"] F [dT(P')«l lnP'] F 


AT(P* ) AT(P' ) <1 In P' d In P" 


terms with higher order 
functional derivatives 
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In practice, it is not necessary to go beyond the second order term, and even 
the first order correction is sufficient for most applications; this will be the 
subject of a paper currently under preparation. 

To summarize, we have called attention to a potentially serious omission in the 
application of linear methods to the retrival of temperatures from satellite 
measured radiances. The weighting functions usually employed do not properly 
take into account that atmospheric transmission itself depends upon temperature . 
We have shown here the form for the proper weighting function. The impact of 
this omission has been shown to be equivalent to a radiance brightness tempera- 
ture „rror of several tenths of a degree for typical tropospheric channels on 
current sounders and as much as 1°K for very narrow channels that may be em- 
ployed in high vertical resolution sounders of the future. More generally, we 
have presented a formalism for correcting transmission functions for changes 
in the temperature profile; the corrections can be computed to any desired 
degree of accuracy. 
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APPENDIX 


We can derive the proper T-weighting function by first dividing the atmosphere 
into isothermal layers and treating the temperature in each layer and at the 
surface as an Independent variable T j , where j = 1 is the topmost layer, j = s - 1 
is the bottom layer, and s corresponds to the surface. The radiance measured 
at the top of the atmosphere can then be written 

s 

ROr T r ... T s ) ■I B(T ) At. (Tj , T 2 , . . . T ) 
j ■ i 

where A t . = r^, - r. , r. is the transmission from the top of the atmosphere 
to the bottom boundary of the jth layer for j < s, r 0 1, and o. The T- 
weighting function is defined as the change in radiance due to a change in the 
temperature in layer i, divided by the temperature change, in the limit as the 
temperature change approaches zero, 



<1B(T) 

dT. 

1 


At. 


* Z] B(T . * 


dAr . 

j 

‘>T, 


where the terms in the summation for j < i are dropped because they are neces- 
sarily zero. 


The above expression is the discrete, or vector form of the T-weighting function. 
The surface component simplifies to 


w s r 


dB(T-) 


dT_ 


s-i 


by virtue of the fact that atmospheric transmission is independent of the surface 
temperature. The atmospheric components can be expressed in continuous form 
by allowing temperature to be a continuous function of a vertical parameter x 
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(ranging from 0 at the surface to oo at the top of the atmosphere). In that case 
the transmission will be a continuous point function of x, and at each point x, it 
will be a continuous functional of T(x') in the range x < x' £ which will be 
denoted by 

t[T(x , )I": x]. 

The radiation at the top of the atmosphere will also be a functional of T(x), 

R[T(x)|£] = B(T s ) t [T(x)|q ; 0] 


f 

•'o 


B(T(x')) t[T(x)|; 

tlx' 


dx' 


The T-weighting function is defined in a manner completely analogous to the 
discrete case, but taking into account the dependence of R on the function T(x). 
We consider a small positive change in temperature AT(x) which is zero every 
where except in the range x - l/2 Ax £ x £ x + l/2 Ax, where A x is a small 
positive quantity. The T-weighting function is then defined as 


W(x) Lim 

AT(x). Ax 


AR [T (x )[ q] 

0 AT(x) | Ax| 


which, by definition, is the functional derivative of R with respect to the function 
T(x), and which will here be denoted by 

d r R IT (x)| q] 

W(x) = -* -• 

[dT (x) dx] F 

(See, for example, Evans, 1964, for a mathematical discussion of functionals.) 
Computing the differentials and then proceeding to the limit, we obtain 


W(x) 


<IB(T) 

dT 


T( * ) 


d_ 

dx 


r T(x' )|"; x] + B(T S ) 


d F r [Tf x ) | q ; 0] 
[dT(x) 1 1 x ] K 


f 

•'o 


(Ldr [T(x)| » ; x ] 

B(T(x')) — - 1 dx' 

[dT(x ) cl x] F • dx 


where the upper limit of the integral has been adjusted in recognition of the 
fact that the integrand is zero for x' > x. 
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We carry this expression into the main part of the paper with the substitution 
dx = -d In P for the height parameter. 
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Figure 1. The T-weighting functions for two of the channels on the NIMBUS 6 
High-resolution Infrared Radiometer Sounder (HIRS). W, the proper T-weight- 
ing function, introduced in this paper, is compared with W ( 1 ) , the T-weighting 
function usually employed in linear methods of retrieving temperature. W REN 
(shown only for Ch. 5) is W renormalized so that its peak value is the same as 
for W". 
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Figure 2. The proper T -weighting function W compared with the usually em- 
ployed weighting function W 11 * for a narrow, practically monochromatic 
channel mid-way between two absorption lines in the CO, 4.3 pm band. 
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